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The small surface tension of the interface between hadronic and quark-gluon-plasma domains, along 
with a negative curvature tension, implies that the uniform plasma is unstable against spontaneous 
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formation of hadronic bubbles. We furthermore show that spherical bubbles are not stable against 
small perturbations. 

1. Introduction 



k>( | In what follows, we shall adopt the common assumption that the quark-gluon- 

plasma phase turns into the hadron phase via a first order transition. In other words, 
that at some temperatures the two phases may coexist. 

First order phase transition usually proceeds by nucleation of bubbles. Consider 
for example, boiling water at atmospheric pressure near the transition temperature 
T = 37 K. The free energy of a bubble of vapor in the water has a free energy which 
depends on its radius: 

F(R) = AP -7iR 3 + a A-kR 2 +... . (1) 

o 

AP is the difference between the pressures of the two phases (zero at the transition 
temperature T = T , and positive above); a is the excess of free energy per unit area 
for a planar interface. 
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F(R) is shown schematically in Fig. 1. For T < Tq, both the volume and surface 
terms act to suppress bubble formation; any bubble formed by a fluctuation collapses 
to R = 0. For T > T , however, the volume term AP47ri? 3 /3 encourages growth 
of bubbles while the surface term a47rR 2 creates a barrier. The maximum of F(R) 
occurs at some R = R c . A bubble formed in the superheated liquid with R < R c will 
shrink away, but a bubble formed with R > R c will grow until the entire system is 
vapor. 
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Figure 1: Schematic representation of bubble free energy in boiling water, retaining only volume 
and area terms, for T < Tg (solid) and T > Tq (dashed). 



In Eq. (1) we have indicated that there might be more terms beyond the area 
term, having in mind an expansion in powers of 1/R. In condensed-matter contexts, 
a "curvature" term, a ■ 8nR, is clearly negligible; by dimensional analysis, a/a ~ A, 
where A is a microscopic scale, and (aR)/(aR 2 ) ~ X/R, an exceedingly small number 
since the typical nucleating bubble is much larger than either the inter-molecular 
spacing or the correlation length which are usually of the order of angstroms. The 
picture in Fig. 1, based on volume and area terms, is then the whole story. 



Note that our picture will be substantially the same for liquid droplets in the 
vapor as for vapor bubbles in the liquid, since a is the same in the two cases. This 
is because the surface tension may be measured in planar interfaces, and it doesn't 
matter which phase is on which side. 

2. The QCD Transition 

In QCD (at least in the limit of massless quarks) there is only one scale. Masses, 

radii, the transition temperature, the latent heat, etc. are all on the order of 200 

MeV or 1 fm. This is then the only "microscopic" size, and 1 fm should also be 

the typical size of a bubble, and one cannot ignore additional terms in (1). Indeed, 

Monte Carlo measurements 2 ' 5 ' 7 and MIT bag model calculations 1 ' 3 found that for 

QCD <j/Tq turns out to be a small number, and one is forced to include additional 

terms, 

4 
F(R) =AP-ttR 3 + a47iR 2 +a87rR + ... . (2) 

o 

Furthermore, it was found 3 that a is sizeable and negative. Fig. 2 depicts F(R) 
for a hadron bubble within the quark gluon plasma. It is qualitatively different than 
that shown in Fig. 1. Indeed, the minimum in the free energy for T > T implies that 
vacuum bubbles are stable. The minimum of F(R) is not a mere artefact of using 
the expansion (2): an exact bag model calculation 3 exhibits a very similar behavior. 

It was further suggested 3 ' 4 that this result might transcend the bag model: All that 
is required is a thin interface, the energy of which is due primarily to the reaction of 
the fields on either side. A lattice calculation of a in the pure glue theory 5 apparently 
supports this picture. 

The fact that one can lower the free energy of the plasma phase by creating a 
vacuum bubble seems at first glance to imply that the equilibrium plasma is a Swiss 
cheese of bubbles, the density of which is determined by the (unknown) bubble- 
bubble interaction. This would only be true, however, if (1) the bubbles do not affect 
each other, and (2) each bubble is stable in itself. A simple calculation will teach us, 
however, that spherical bubbles are unstable against deformation, and that growth 
of long, narrow structures is favored. We do not address the problem of determining 
the actual equilibrium configuration of the high-temperature phase, a problem which 
will require much more sophisticated analysis. 

In order to study the stability of the spherical bubble against small perturbations, 
we expand the free energy around the equilibrium spherical bubble, of radius Ro- 
Minimizing (2), we find that Rq satisfies 

4ttAP Rl + 8naR + 8na = 0. (3) 

The generalization of (2) to non-spherical shapes is 6 

F = APV + aJdS + ajdsQ- + ^-)+---, (4) 
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Figure 2: Free energy of a hadron gas bubble in the plasma for temperatures near 
the transition temperature, fixed at T = 150 MeV. T = 155 MeV (solid), T = T = 
150 MeV (dashed) and T = 147 MeV (dotted). The figures are taken from a bag 
model calculation of Ref. 3. 



where the integrand in the third term is the local extrinsic curvature, with R\ and 
R 2 representing the principal radii of curvature. Let the bubble surface be specified 
by 

R(9,0)=R O + 5R(9,0). (5) 

The total 0(SR) contribution to F vanishes identically on the equilibrium radius 
R = Rq. The correction to the free energy is then second order and may be written 
as 

5F= I sin 6 d9 d(p 5R (A + BL 2 ) 5R. (6) 

— L 2 is just the angular part of the Laplacian operator, and 

A = APR Q + a, B = ?- + ^-. (7) 

Recall that we presume a < 0. When R is determined by (3), it is easy to see that 
A is positive and B negative. 

Expanding 5R in a multipole expansion, 

<$R(M) = E c r%»(M), (s) 

l.rn 

we have 

5F = Y J \A + Bl{l+l)]\cT\ 2 . (9) 

l.rn 

Evidently SF > for a monopole perturbation, since (2) has been minimized; for a 
dipole perturbation we obtain SF = , consistent with the fact that I = 1 corresponds 
to a translation of the bubble. For quadrupole and higher deformations, however, SF 
is always negative, showing instability against arbitrary changes of shape. 

3. Discussion 

We have here considered only static properties of the quark-hadron interface. 
Dynamical processes in the course of the phase transition will destabilize the interface 
even more. For example, the necessity of carrying heat and baryon number away from 
a growing vacuum bubble will lead to a negative effective surface tension and hence to 
complex phenomena of pattern formation. On the other hand, interaction of bubbles, 
which was completely ignored in all the above treatments, may play a crucial role in 
stabilizing the bubbles' surface. 

A word of caution: the above analysis relies on a bag model calculation 3 . And 
although Monte Carlo measurements confirmed the smalleness of a as well as the 
sign and magnitude of a of the bag model calculation, it is still far from actually 
confirming even qualitatively the appearance of a minimum in F(R) for a hadron 
bubble with finite radius surrounded by the quark gluon plasma, as depicted in fig. 
2. As a matter of fact, the lattice results show that the interface is about 2fm thick, 



and it is clearly not negligible compared to the temperature or to Rq , whereas the 
stability analysis is based on a thin wall approximation. 
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